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ABSTRACT 

We present semi-analytical similarity solutions for the inside-out, expansion-wave 
collapse of initially virialized gas clouds with non-isothermal equations of state. 
Results are given for the family of negative- index polytropes (P oc p 7 , 7 < 1), but 
we focus especially on the so-called logotrope, P/P c = 1 + A ln(p/p c ). In a separate 
paper, we have shown this to be the best available phenomenological description of 
the internal structure and average properties of molecular clouds and dense clumps 
of both high and low mass. The formalism and interpretation of the present theory 
are extensions of those in Shu's (1977) standard model for accretion in self-gravitating 
isothermal spheres: a collapse front moves outwards into a cloud at rest, and the gas 
behind it falls back to a collapsed core, or protostar. The infalling material eventually 
enters free-fall, so that the density profiles and velocity fields have the same shape 
(p oc r~ 3//2 and — v oc r~ 1//2 ) at small radii in logotropic and isothermal spheres both. 
However, several differences arise from the introduction of a new equation of state. 
The accretion rate onto a protostar is not constant in a logotrope, but grows as 
M oc t s during the expansion wave. Thus, the formation time for a star of mass M 
scales as M 1 / 4 ; low-mass stars are accreted over longer times, and high-mass stars over 
shorter times, than expected in isothermal clouds. This result has implications for the 
form and origin of the stellar IMF. In addition, the gas density behind an expansion 
wave increases with time in our theory, but would decrease in an isothermal sphere. 
The infall velocities also grow, but at an initially much slower rate than found in an 
isothermal collapse. These results apply to low- and high-mass star formation alike. 
We briefly discuss how they lead to older inferred collapse ages for Class protostars 
in general, and for the Bok globule B335 in particular. 

To appear in The Astrophysical Journal 



1. Introduction 



One of the cornerstones in any theory of star formation has to be a correct description of 
the gravitational collapse of gaseous clouds. The pioneering studies of Bodenheimer & Sweigart 
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(1968), Larson (1969), and Penston (1969) gave numerical integrations of the isothermal collapse 
of uniform-density spheres, and Larson and Penston both went so far as to derive semi-analytical 
similarity solutions (descriptions of a fluid flow in which the densities and velocities approach 
invariant forms) for the evolution of their spheres. In a seminal paper which defined the current 
paradigm for "inside-out" collapse and low-mass star formation, Shu (1977) derived a different, 
one-parameter family of similarity solutions for the isothermal problem. Hunter (1977) found still 
more, and Whitworth & Summers (1985) finally pointed out that there is in fact a two-dimensional 
continuum of such solutions. Common to all of these is the eventual formation of a central point 
mass which subsequently grows as the cloud continues to fall onto it. This core is naturally 
identified as a protostar (or star and disk). 

One of Shu's (1977) solutions in particular — the so-called expansion wave — has played a 
central role in the development of the standard theory of star formation, as reviewed, for example, 
by Shu, Adams, & Lizano (1987). In part, this is because of the simple and intuitive nature of the 
expansion wave, and because it lends itself well to the (perturbative) inclusion of rotation (Tereby, 
Shu, & Cassen 1984) and mean magnetic fields (Galli & Shu 1993a, b) in isothermal clouds. In 
addition, however, it is physically relevant because it describes the relatively gentle collapse, from 
the center on out, of a sphere which is initially in virial equilibrium. As we discuss shortly, the 
idea that pre-collapse interstellar clouds are in equilibrium configurations is one which meshes 
with a variety of observations. There are contrary arguments (e.g., Bonnell, Bate, & Price 1996) 
that star-forming clouds should collapse from non-equilibrium, or even near-uniform conditions; 
but this is a violent process (the flow of Larson [1969] and Penston [1969] is everywhere highly 
supersonic) which would give rise to strong kinematic signatures that have never been observed. 
Rather, the most convincing claim to a direct detection of collapse has been made for the Bok 
globule B335, where molecular line profiles are consistent with the more moderate velocity field 
expected for an expansion wave (Zhou 1992; Zhou et al. 1993). 

The standard theory has its difficulties, however. For instance, Shu's (1977) model describes 
the collapse of highly unstable, singular isothermal spheres, and is not necessarily applicable to 
more realistic clouds with finite central densities (such as the marginally stable Bonnor-Ebert 
sphere: Hunter 1977; Foster & Chevalier 1993). As another example, an isothermal expansion 
wave leads to a time-invariant protostellar accretion rate, M = 0.975cr 3 /G, which appears to 
be inconsistent with the observed luminosities of young stellar objects in the Taurus-Auriga 
star- forming region specifically (e.g., Kenyon et al. 1994). But most important is the fact that 
an isothermal equation of state, while a natural and instructive first approximation, is ultimately 
incompatible with observations of both giant molecular clouds (GMCs) and the subcondensations, 
the dense clumps where stars are ultimately born, within them. 

As suggested above, GMC complexes, and individual clumps, are well described by models of 
(magnetic) virial-equilibrium spheres which are "truncated" by a constant surface pressure exerted 
by a diffuse external medium (e.g., Myers & Goodman 1988; Bertoldi & McKee 1992; Elmegreen 
1989; McLaughlin h Pudritz 1996). (The dense clumps are often referred to as cores, but here we 
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reserve that term for the central point mass, or protostar, which develops in a collapsing sphere.) 
The virial balance relies on superthermal random velocities which increase outwards within any 
given clump. The nonthermal component ctnt °f the total velocity dispersion (linewidth) is larger 
in higher- mass clumps (Caselli & Myers 1995), but is present at even the smallest observed radii in 
low-mass clumps as well (Fuller Sz Myers 1992). Indeed, (tnt grows more rapidly with increasing 
radius in low-mass clumps than in high-mass ones, and no interstellar cloud obeys an isothermal 
equation of state. 

The first attempts to model the nonthermal linewidths (e.g., Myers & Fuller 1992; Lizano & 
Shu 1989) tended to focus on the fact that ctnt is small in low-mass clumps. Myers & Fuller, 
in particular, construct equilibrium, self-gravitating spheres from a gas equation of state (EOS) 
which is explicitly isothermal at small radii in a cloud, and becomes progressively "softer" (P 
depends less than linearly on p, so that cr 2 = P/p increases with decreasing density) further out. 
Because the clumps in GMCs can be viewed as pressure-truncated spheroids, with truncation 
at a smaller radius corresponding to a lower mass, this approach effectively ensures that the 
low-mass clumps are still treated as (very nearly) isothermal spheres. Any consequences which 
the nonthermal support might have for collapse and star formation are then largely confined to 
the high-mass regime. For example, Myers & Fuller (1992) estimate the formation times of stars 
in their model, and conclude that high-mass stars are accreted relatively faster than they would 
be under an isothermal EOS, while low-mass stars grow at roughly the same rate. 

By contrast, McLaughlin & Pudritz (1996=MP96) have adopted a phenomenological, 
logotropic EOS — P/P c = 1 + A ln(p/p c ) — for molecular clouds and clumps. This EOS provides 
a unified treatment of low- and high-mass clumps and entire GMCs. It can account quantitatively 
for (1) the global size-linewidth and mass-radius relations (e.g., Larson 1981) between cloud 
complexes; (2) the typical mass and density contrasts between GMCs and the clumps within 
them; and (3) the observed dependence of linewidth on radius inside dense clumps of any mass. 
Moreover, the form of the EOS itself is never close, even on small spatial scales or in low-mass 
clouds, to the isothermal P oc p. The details of both low- and high-mass star formation in a 
logotrope will therefore differ from the standard results for isothermal collapse. (It should be noted 
that our logotrope is different from that introduced by Lizano & Shu [1989], which essentially has 
P ~ p + hip, and cannot account for all three of the properties of interstellar clouds listed here.) 

In this paper, then, we derive semi-analytical similarity solutions for the collapse of 
nonmagnetic, non-rotating gas spheres with softer-than-isothermal equations of state. We 
concentrate on the case of a logotropic EOS as the most self-consistent description of the 
ubiquitous turbulence in interstellar clouds; but in a pair of Appendices, results are also given for 
the more general family of negative-index polytropes (P oc p 7 , with < 7 < 1). (Although the 
self-similar collapse of polytropes has previously been investigated by, e.g., Suto & Silk [1988], we 
also briefly discuss the accretion rates and star formation timescales in such spheres.) 

Section 2 begins with a brief summary of the equilibrium structure of a logotrope (MP96), 
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and then goes on to an analysis of its collapse. Our development owes much to the physical 
discussion of Shu (1977): we first find a class of collapse solutions for singular, p oc r _1 spheres 
which are very accurate approximations to more realistic, non-singular logotropes, and then focus 
on the inside-out, expansion- wave collapse of an initially virialized cloud. Two main features are 
a post-collapse density profile and velocity field which have the same free-fall shape (p oc r -3 / 2 
and —u oc r -1 / 2 at small radii) found for an isothermal expansion wave, but a rate of accretion 
onto the collapsed core which is no longer constant in time; rather, M oc i 3 here. These results 
allow for an evaluation, in §3, of the formation times for stars of various masses. High-mass stars 
will be accreted more rapidly in a logotrope than in an isothermal sphere, but a low-mass star in 
a logotrope will appear later than a star of the same mass would in an isothermal sphere. This 
"squeezing" of formation timescales, which we estimate to be typically of order 1 — 3 x 10 6 yr, 
should have important implications for the origin and shape of the stellar IMF. In §4 we consider 
the collapse of a 1 M logotrope specifically, and briefly discuss relevant observations of two 
young stellar objects, NGC 2071 and B335. We conclude that the age of B335 is probably greater 
than 10 6 yr, an order of magnitude older than suggested by the application of isothermal collapse 
models. Since B335 is one of the family of the youngest known, "Class 0" protostars (Andre, 
Ward-Thompson, & Barsony 1993), these objects might generally be rather older than is currently 
believed. 



2. Similarity Solutions 
2.1. The Logotrope 



P l + Alnf-^-) , (2.1) 



Our logotropic EOS (MP96) is 



Pc \Pc 

where P c and p c are the central pressure and density in a cloud, and A is essentially a free 
parameter. The total velocity dispersion a 2 = P/p increases outwards from r = until it reaches 
a maximum at some large radius, beyond which it decreases somewhat. Since linewidths are 
essentially purely thermal on the smallest scales in real interstellar clouds, and are also observed 
to increase with radius (e.g., Fuller & Myers 1992; Caselli & Myers 1995), the central velocity 
dispersion in a logotrope is identified with the thermal value: a 2 = P c j p c = kT/pmu- Once a 



value of A is specified, the EOS (2.1) can be used to integrate a Lane-Emden type of equation 
for a self-gravitating, virialized gas sphere. This then allows for a quantitative comparison of the 
internal a vs. r profiles of real and model clouds. Good agreement with data on GMC clumps, of 
low and high mass alike, is obtained for A ~ 0.2 ± 0.02. 

These points are discussed in more detail in MP96. Starting from the premise that interstellar 
clouds are pressure-truncated — their boundaries defined by a surface pressure P s due to a more 
diffuse, surrounding medium — that paper develops expressions for the masses and densities of 
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virial-equilibrium gas spheres with arbitrary equations of state. Quite general arguments then 
imply that GMCs, and the largest clumps within them, are all at or near their critical masses: 
if they were any more massive, given their temperatures, surface pressures, and magnetic field 
strengths, they would be unstable to radial perturbations (see also McKee 1989; Bertoldi & McKee 
1992). Although this conclusion is independent of the exact form of the gas EOS, not all equations 
of state provide for such a critical mass. For example, the family of negative-index polytropes (in 
which P oc p 7 , with 7 < 1) also have velocity dispersions that increase with radius, but are always 
stable against radial perturbations in the analysis of MP96. Put another way, and independently 
of any stability analysis, such polytropes cannot simultaneously account for the internal structure 
of molecular clumps and the global properties (Larson's [1981] laws: M oc R 2 and <r ave oc R 1 / 2 ) of 
GMCs; the logotrope is the simplest EOS for which this is possible. 

The critical mass of an A = 0.2 logotrope is given by equation (2.11) of MP96: 

M -it = J| M (j^j ( L3xl0 5^ cm -3 K ) ' (2 ' 2) 

where we have explicitly used the fact that a 2 = kT/ixmu, and set the mean molecular weight 
to ji = 2.33. The constant a here depends on the strength of any mean magnetic field threading 
the cloud. If there is equipartition between kinetic and mean magnetic energies, then a ~ 1; in 
the absence of any mean field at all, a = 35/18 for the logotrope and M cr j t ~ 92M for this 
temperature and surface pressure. A smaller value of A results in a larger critical mass: for 
example, A = 0.18 would raise M cr it by a factor of about 2.5. MP96 also derive the mean density 
of a critical-mass logotrope: 

" 131 * ltr " g Cm " (line)" ( l.3xl/tcm-»K ) ' (2 ' 3) 

The fiducial P s chosen here is appropriate for dense clumps embedded in self-gravitating 
GMCs, which amplify the pressure of the hot ISM by about an order of magnitude (Elmegreen 
1989; Bertoldi & McKee 1992); it is also similar to the value adopted by Shu (1977) in his study 
of isothermal collapse. It is significant that the logotropic M cr jt is vastly greater than the ~1Mq 
mass of a critical, T = 10 K isothermal sphere under the same surface pressure. This is part of 
what makes our model more consistent with observations of the ISM, and in particular with the 
arguments, mentioned above, that only the largest clumps in GMCs are at their critical mass, 
with smaller ones being in a stable virial equilibrium. 

Because we are dealing with pressure-truncated spheres, a cloud's mass can be viewed as 
being set by the dimensionless "truncation radius" R/r (where r = Sac/lAirGpc^ 1 ! 2 ) at which the 
internal pressure equals that of an ambient medium. That is, there is an EOS-dependent R/ro 
which corresponds to a critical-mass cloud; any smaller R/ro is equivalent to a unique M/M cr it < 1. 
Indeed, R/ro, along with a surface pressure P s and kinetic temperature T, determines all the 
properties of a cloud (see Appendix A of MP96). An example of this is given in Table [l], where we 
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list the center-to-edge density and pressure contrasts, mean density, and mass-averaged linewidth 
(thermal plus nonthermal: cr 2 ve = P ave /Pave) for logotropes with A = 0.2. The absolute values of 
the masses and free-fall times there (tg = [37r/32G / o ave ] 1 / 2 ) apply to purely hydrostatic spheres 
with no mean magnetic fields, as these are the types of objects which our collapse calculations will 
directly address. If P s or T is changed, the mass and mean density of the critical cloud will be 
affected as in equations (|2.2| ) and fl2.3|) . Then, because the dimensionless R/ro actually determines 
the mass and free-fall time relative to the critical values, these quantities will both change at every 
radius in Table [j], and by the same factors as in the critical cloud. In any case, it is significant that 
the total linewidth a ave is generally larger in bigger clouds. Lower-mass objects are supported to 
a greater, but not complete extent by thermal motions. 

Figure [l] shows the density profile of a critical-mass logotrope, and compares it with that of a 
critical-mass, or Bonnor-Ebert, isothermal sphere (Bonnor 1956; Ebert 1955; McCrea 1957). The 
structure of less massive spheres can also be derived from this diagram; for example, the p vs. r 
curve for a logotrope with M/M cr ; t — 0.33 is obtained by vertically shifting the curve in Fig. [I] 
such that pj p s = 1 at r/ro ~ 9.4 (cf. Table |l]). Also shown in the Figure, as broken lines, are 
the analytical singular solution to the equation of hydrostatic equilibrium for a self-gravitating 
logotrope (MP96), 

p(r) = {AP c /2nG) 1 ' 2 r- 1 , (2.4) 

and the singular solution p(r) = (a 2 /27rG) r~ 2 for an isothermal sphere. Equation ([O]) closely 
approximates the bounded logotrope (finite central density; solid line) at all r/ro ^ 0.4, and 
thus adequately describes the structure even of stable low-mass spheres. On the other hand, the 
singular isothermal sphere bears very little relation to the Bonnor-Ebert configuration. (The 
logotropic prediction p oc r _1 is one with some observational support. This is discussed in MP96; 
here we note only that an observed column density profile TV oc r , which is not uncommon, does 
not necessarily imply p oc r -2 .) 

Since a low-mass clump (say, M tot ~ 1M ) is expected to be globally stable against radial 
perturbations, how does it come to collapse and form a star? When magnetic fields are present, 
ambipolar diffusion should induce a slow contraction of the clump, during which its initial (time 
t = — oo) central concentration will be enhanced and its density profile brought into agreement 
with the singular solution ( |2.4| ) over essentially all radii — a configuration that is manifestly 
unstable (e.g., MP96). Alternatively, in a purely hydrodynamic situation, collapse could be 
initiated by the growth of a local density fluctuation; or the clump could be put out of equilibrium 
by the accretion of mass from the surrounding medium. The time at which a cloud becomes 
singular is labelled t = 0, and corresponds to the formation of a central point mass (or core; 
physically, a protostar and disk). Following Shu (1977), it is the subsequent, t > evolution of 
nonmagnetic spheres that we now examine in detail. 
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2.2. Fluid Equations 



We begin with the usual Eulerian equations for spherically symmetric flow in a self-gravitating 
fluid (cf. Larson 1969; Shu 1977): 



du du GM dPdlnp 

dt dr r 2 dp dr 

dM A 2 
— - + Airr pu 
at 

4irr p 



Or 






. 



(2.5) 
(2.6) 
(2.7) 



where M(r, t) is the total mass (including any core) inside radius r at time t, and u(r, t) is the 
instantaneous fluid velocity. The gas EOS influences the collapse through the pressure force 
p~ 1 dP/dr in equation (2.5). For the most part, we confine our attention to the logotropic 
EOS (eq. p.lj ), but in Appendices A and B below we also consider the more general family of 
negative-index polytropes. 



We define a similarity variable 



r 

a t t ' 



(2i 



where at is a function of time with dimensions of velocity. Dimensionless density, velocity, and 
mass variables may then be set up as 

a(x) = 4irGt 2 p(r,t) v(x) = u{r,t)/a t m(x) = GM(r,t)/<4t . (2.9) 

Further, because the square of the sound speed is dP/dp = AP c /p for this EOS, we are led to set 

a t = [AP^Gt 2 )} 1 / 2 , (2.10) 

where A ~ 0.2 for real interstellar clouds, and where the constant P c refers to the central pressure 
in the initial, pre-collapse (t = — oo) cloud (as such, it reflects the constant pressure of the ambient 
medium, and the total mass of the cloud; see, e.g., Table p]). Ultimately, the central core mass 
grows as M(0,t) oc a^t oc t 4 , and the accretion rate as M(0,t) oc i 3 . This is in sharp contrast 
to an isothermal collapse, where the accretion rate is constant in time (see §3 below for further 
discussion). In any case, straightforward manipulation of equations (|2.5|) - (2.7) yields 



(2x - vf 



d lna 



dx 



x 



-(2a: -v) 



(2x -v) + [2(2x -v)+v] 



(2x -v) 2 

a 



dv 
dx 

m 



a , ,2 

(2x-v) 

4 ax 



(2x -v) + 



- + v(2x 
a 



ax 



(2x - v) 



(2.11) 
(2.12) 
(2.13) 



These equations are written in a form similar to that of Shu (1977) for the isothermal sphere, 
which is included as a special case in the general development of Appendix A below. 
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There are two simple analytical solutions of equations ( |2.11 ) - (2.13). The first has a uniform 
density, 

2 2x 2x 3 . A , 
a = - v = — m = , (2.14) 

3 3 9 v ' 

and is equally valid for any polytropic EOS. Of more interest is the static (singular) solution, 

V2 x 2 
a = — v = m=—^. (2.15) 
x ^2 

This corresponds to a time-independent density profile, p(r) = (AP c /2ttG) 1 ^ 2 r" 1 , which is just 
that found (eq. |2.4| ) to satisfy the equation of hydrostatic equilibrium for a self-gravitating 
logotrope. As discussed in §2.1, it also describes the large-radius structure of the bounded sphere 
— the type of cloud we expect to begin with at t = — oo. 

An entire class of singular solutions, valid at large radii or small (positive) times, can be 
found by taking a and \v\ to be much smaller than x in equations ( |2.11| ) - ( |2.13| ). In this x — > oo 
limit, we have 

C (7/ 2\ Cx 2 



„^__^__j m ^_. (2 . 16) 

For an inflow scenario (v < 0), the density normalization C must be greater than the hydrostatic- 
equilibrium value \[2. Each C > \[2 therefore corresponds to the collapse of a logotrope which is 
initially just out of equilibrium, with densities (and pressures) everywhere enhanced by a factor 
1 + 5 = C j\[2. As in the equilibrium case, the cloud starts at t = — oo with a finite central density, 
and p at large radii given by equation ( p.!6| ), and then contracts to follow that singular profile 
everywhere by t = 0. These solutions are the analogues of Shu's (1977) "minus solutions without 
critical points" for the singular isothermal sphere — although we show in §2.3 that in our case, 
some of the out-of-equilibrium flows can pass through critical points. 

Although v is a finite constant at t = here, the dimensional velocity u ~ vt vanishes, and 
the cloud is essentially at rest at the time of core formation. This is consistent with a physical 
picture at t < of quasistatic evolution towards the singular profile. As suggested above, such 
a slow contraction could be driven by ambipolar diffusion, and we therefore view the family of 
solutions (|2.16| ) as reasonable starting points for the post-core formation evolution of a collapsing 
logotrope. This is the same tack taken by Shu (1977) in his treatment of the isothermal sphere, 
but in one important way it is perhaps better justified for our EOS. If it happens that the mean 
magnetic field is initially weak and unable to significantly influence the evolution of a cloud, then 



the assumption of a gentle approach to the singular profile, and the validity of equations (2.16) in 
the t = limit of the collapse, have to be questioned 

The singular isothermal sphere is a poor description of the critical-mass Bonnor-Ebert 
configuration, or of any less massive one (see Fig. |l]). Thus, purely hydrodynamic collapse 
simulations of the latter clouds (e.g., Hunter 1977; Foster & Chevalier 1993) differ significantly 
from the results of Shu (1977). The approach to an r~ 2 profile at t = is rather violent for them, 
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and involves appreciable fluid velocities that invalidate the main premise (u = at t = 0) of the 
similarity solutions discussed here and by Shu. The subsequent t > collapse then proceeds much 
more rapidly than for a singular sphere, with an accretion rate that decreases over time. In our 
case, the logotropic singular solution (p oc r~ x : eq. | 2.4 ] and Fig. [j]) does describe the structure of 



bounded clouds of any interesting mass. The adjustment to a fully singular profile at t = should 
then be a relatively easy one, and it is reasonable to expect that equations ( 2.1 6|) adequately 



describe any collapsing logotrope at the point of core formation. This argument is consistent with 
numerical simulations (Foster & Chevalier 1993) for the collapse of a bounded isothermal sphere 
whose central concentration is ~70 times that of the Bonnor-Ebert value. Although such a cloud 
is highly unstable, it is well approximated at large radii by the singular isothermal sphere, and its 
accretion rate quickly approaches the constant value predicted by Shu (1977). 

Turning now to the opposite limit of small radii, x —* 0, we have that \v\, a — > oo for collapse, 
and thus 

a -»■ 4 — -r v -»■ - m^m , (2.17) 



\2x 3 J 



.r 



where the reduced mass mo of the collapsed core also serves to fix the accretion rate onto the core 
(Mo oc mot 3 by eq. [ |2.9| ; see also §3 below). The effects of shocks (e.g., Tsai & Hsu 1995) and 
radiation pressure (for massive stars; Jijina & Adams 1996) will invalidate this solution at very 
small x, but should not affect our discussion of the area external to such disturbed regions, nor 
the value of mo itself. 

The x~ 3 / 2 and x _1 / 2 scalings of a and v here are identical to the corresponding results for the 
isothermal sphere (Shu 1977), or any other polytrope with < 7 < 1 (Appendix A). This reflects 
the fact that the ratio of gravitational and pressure forces, 



GM/r ma 



p 1 dP/dr x \dlnxj 3 x 5 / 2 



1 A./o ™ 3 / 2 



m 







(2.18) 



increases without limit as x — > 0, and the gas at very small radii approaches the core in free-fall. 
At any fixed time t > 0, the dimensional velocity field near the center therefore tends to 
u = (2GM0/V) 1 / 2 , and the density profile to 

M _ M _ 3//2 
P{r} = 4^ = 47r(2GM )V2 r ' (2 ' 19) 

which is equivalent to equation ( |2.17| ) for a spatially invariant accretion rate M. When written 
this way, the effects of the gas EOS are completely contained in the values of M and Mo at a 
given time; at least in the absence of rotation and magnetic fields, observations of the shape of 
the density profiles or velocity fields in the central regions of a collapsing cloud are insufficient 
to constrain the underlying equation of state. Any data which favor an r~ 3 / 2 density profile in 
such objects (e.g., Butner et al. 1990, 1991) are qualitatively consistent with both the logotropic 
hypothesis and the isothermal one (quantitatively, however, the different accretion histories in 
the two models may well allow for an empirical discrimination between them; see §4 below). By 
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the same token, observations which suggest a different structure (e.g., Andre et al. 1993) are a 
challenge to both scenarios. 

A full similarity solution for the collapse of overdense clouds, over all times t > and all 
radii r > 0, begins at large x as in equation ( |2.16| ), and eventually tends to the form ( 2.17| ) at 



small x. As such, there is a one-to-one relation between C and uiq. This is obtained by numerical 
integration of equations (2.11) and ( |2.12j ), and is laid out in Table 0. Evidently, the larger the 



overdensity C/y/2, the larger is the reduced core mass mo, and the faster the accretion rate. The 
limit C — ► \pl refers to the collapse of a logotrope which is initially arbitrarily close to hydrostatic 
equilibrium. In this situation — which is the parallel of Shu's (1977) expansion-wave solution for 
the singular isothermal sphere — a collapse front propagates outwards in a cloud at rest, with the 
material behind it falling away to the center. 



2.3. Critical Points and the Expansion Wave 



A look at equations ( [2.11| ) and ( 2.12| ) shows that when a collapse solution passes through a 
point x* > which has 

2x* - v* - o~ 1/2 = , (2.20) 

the derivatives da/dx and dv/dx will be undefined unless the quantities on the right-hand side of 
the fluid equations also vanish there. Any valid flow must therefore satisfy 



1/2 . 

aj 1 



+ + 2x^ = (2.21) 



at such critical points. Mathematically, the a(x) and v{x) which everywhere satisfy these two 
relations would present a third analytical similarity solution for our problem; but it is one which 
relies on a rather special interplay between the fluid variables, and one which we therefore neglect. 
Instead, physical flows must either avoid the line of critical points (or singular points, or sonic 
points) defined by equations ( |2.20| ) and ( 2.21 ), or else pass right through it. 



Whitworth & Summers (1985) elaborate on the topology of critical points and the implications 
for a self-similar, isothermal collapse, and we refer the interested reader to that paper for an 
excellent discussion of such issues. Here we note that, for any gas EOS, a linearization of the fluid 
equations in the neighborhood of any critical point x* leads to three pairs of eigenvalues for the 
slopes da/dx and dv/dx there, and hence to three eigensolutions for a(x) and v(x). One of these 
is physically inadmissible but the other two, corresponding to the "minus" and "plus" solutions of 
Shu (1977), have the velocity either decreasing or increasing towards smaller x. They are described 
quantitatively, for the logotrope and for negative-index polytropes, in Appendix B. 

It is the minus solutions, for which v — > — oo as x — * 0, that are of most interest in the collapse 
context. The velocity fields for some representative examples are shown as the solid lines in Fig. [2|. 
Also shown there are the locus of critical points (dash-dot line) defined above, and a few plus 
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solutions (dashed lines). We will not discuss the plus solutions, other than to say that (1) their 
time-reversed counterparts could be of interest in a wind problem, and (2) in the limit x — ► oo, 
they have a ~ x -4 , v ~ x 2 , and av 2 = pu 2 /AP c = 2. 

Our Fig. Q should be compared to that of Shu (1977). The heavy black line here is the 
expansion-wave solution mentioned above. It passes through the critical point x* = 0.02439, and 
is the limit of two different sequences of minus solutions. Those which pass through critical points 
x* < 0.02439, to the left of the expansion wave, go on to a stagnation point (v = 0) at some finite 
x, beyond which there is outflow. These solutions are therefore of no interest to us, although they 
might be useful in modelling the effects of accretion shocks (see Tsai &: Hsu 1995). If they were to 
be simply truncated at the point where v = 0, they would have vanishingly small spatial extent at 
t = 0, which is unacceptable; and, as in the isothermal case, the static solution (eq. p,15 |) cannot 



be used to continue such flows to infinite x because there is a mismatch in the densities at the 
stagnation point. The densities and velocities still follow the power laws of equation ( |2.17| ) at 
small x, though, with a core mass mo that is uniquely defined by the location of the critical point 
x* and is always less than the expansion-wave value. 

Fundamentally, the expansion wave is identified as an extreme member of this family of 
solutions because its stagnation point is also a critical point, x* = x ew . In this case (see Appendix 
B), 

1 1 

a ew = 8 v ew = m ew = y= mo = 6.67 x 10 . (2.22) 



ew 4^2 6W 6W CW 32^2 

These results do agree with the static cloud profile (|2.15| ) at x = x ew , and the core mass mo 
corresponds to that in the C — > y/2 collapse (cf. Table || above), so an extension to arbitrarily 
large x is obtained by simply joining the two solutions in the obvious way. In fact, the logotropic 
expansion wave comes into x ew from the left with velocity and density derivatives, (dv/dx) Xcvi = 
and (da/dx) Xcv = — 1/32\/2, which also agree with those of the static profile, i.e., the flow is 
continuous at x ew - This is the case for any polytropic EOS with 7 < 3/5, but when 7 > 3/5 
(including the isothermal 7=1) there is a discontinuity at the matching point (see Appendix B). 

Physically, as discussed above and by Shu (1977), an expansion wave corresponds to the 
inside-out collapse of a cloud which is in hydrostatic equilibrium at the moment t = of core 
formation. At any point in time (before the expansion wave reaches the initial cloud surface), a 
total reduced mass of m ew is involved in the collapse, with a small fraction mo/m ew — 0.03 of that 
having already fallen onto the central core. It may also be shown (for any gas EOS) that since the 
point x ew is fixed in similarity space, the wave propagates into the cloud always at the local speed 
of sound, {dP/dp) 1 ' 2 . 

However, the logotrope differs from the isothermal sphere, in that its expansion wave 
passes through two critical points: x* = 0.02439 and x* = x ew . (For an isothermal EOS, the 
expansion- wave v is just tangent to the locus of critical points at x eW) and does not cross it 
anywhere.) It is the presence of these two distinct critical points, both with v* < 0, which also 
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marks the expansion-wave flow as the limit of a second sequence of minus solutions, namely, the 
out-of-equilibrium collapse profiles encountered in §2.1 (eq. [2.16|). Several of these are drawn 
on Fig. ||, to the right of the expansion wave, and labelled by their density normalizations C. 
Evidently, for C close enough to the hydrostatic-equilibrium value of v2, such solutions — which 
have v < everywhere — also pass through two critical points. This is again in contrast to 
the isothermal situation, in which all out-of-equilibrium collapses are free of critical points. On 
some level, however, this is just a mathematical issue; there is nothing intrinsically unphysical or 
artificial about flow through a critical point.Q 



3. Accretion Timescales 

Once numerical solutions to the self-similar fluid equations have been obtained, and the 
reduced mass m,Q found for any minus solution, we have a description of the growth of the central 
core in a collapsing cloud. The formation times for stars of various masses can then be computed. 
Again, we focus on the results for logotropic expansion waves in this Section. Generalizations 
to out-of-equilibrium collapses, and to more general polytropic equations of state, are fairly 
straightforward (cf. Appendix A). 

From equations ( p.!3| ) for the dimensionless mass m and fl2.17| ) for our similarity solutions at 
small radii, the central core mass is given by 

M(0,i) = M 



Here p c and p ave are the central and mean densities of the initial (t = — oo) cloud; tg is its mean 
free-fall time, (3ir/32Gp avc ) 1 /' 2 ; an d it will be recalled that A = 0.2 in the logotropic EOS is most 
appropriate for real molecular clumps and GMCs. The ratio p c /pa,ve and the free-fall time are 
referred to their values in the non-singular cloud because these are then related to the surface 
pressure, kinetic temperature, and total mass (by, e.g., Table [l| above), which fundamentally 
control the collapse. As also discussed in §2.1, a c is just the thermal linewidth, (kT/ /imjj) 1 / 2 . 



1 Fig. |2j shows that the C ~ 1.4245 collapse just brushes the locus of critical points at a single x» ~ 0.1034. 
Solutions for C > 1.4245 then bypass all critical points, while those with v2 < C < 1.4245 encounter two critical 
points (one less than 0.1034, and one greater). In the language of Whitworth & Summers (1985), the points with 
x, < 0.1034 are saddles; those with x* > 0.1034, nodes; and x, = 0.1034 itself is degenerate. (The direction of 
stable numerical integration is away from a saddle, but towards a node; and only a pure plus or minus solution can 
pass through a saddle, while any function can cross a node.) The degeneracy moves towards larger x for increasing 
polytropic index 7, and occurs at x* = 1 = x cv , in the limit 7 = 1. This is why the isothermal expansion wave is 
only ever tangent to the line of critical points, and why the overdense collapse solutions (C > 2 in that case) do not 
involve critical points at all. 



[AP C (4*G)\ 
'3tt 2 



3/2 m o^ 



A 3/2 



G 

p c \ r mod* 



Pa 



t 



if 



G 



(3.1) 



-13- 



The logotropic EOS already accounts for the presence of turbulence in the initial cloud, and a c is 
therefore not to be interpreted as the type of "effective" sound speed discussed by, e.g., Shu et 
al. 1987. The accretion rate is therefore 

where (Mo) is the mean rate Mo/t. As we discuss further in §4, the time dependence of this 
quantity is the most important consequence of the move away from an isothermal EOS, which 
gives a constant Mo = (Mo) = 0.975crj?/G for the expansion wave (Shu 1977; also Appendix A 
below). (Mo) increases with time in our case, and in any negative-index polytrope, because density 
profiles shallower than r -2 put more mass at larger radii in a cloud. 

Now consider a pressure-truncated logotrope which is in hydrostatic equilibrium, and large 
enough that the singular r" 1 density profile holds at the edge of the cloud. The analysis of MP96 
gives the total mass and radius as 

M ^^(^\d (3.3) 

\ Pave J Cr 

and 

2(3A)V2 / Pr \V2 



R = v ' — iff a c . (3.4) 

7T VPave/ 

Apart from corrections for the presence of mean magnetic fields (MP96), these expressions are 
valid for a cloud of any mass as long as the r _1 profile is achieved inside it. In turn, as we have 
seen (cf. Fig. [l]), this occurs for the full range of clump masses — from < 1 to several hundred 
Mq — seen in GMCs. As the expansion wave moves outwards in a logotrope, equations fl3.1| ) and 
give for the fractional core mass (independently of the EOS parameter A) 

[=-) mo, (3.5) 



M to t 16 V2 W 

where, again, too = 6.67 x 10 -4 . Also, the definition of at (eq. |2.10| ) locates the boundary of the 
cloud in similarity space: 

x=*=*4{iy. (3.6) 

a t t 7T 2 \tgj ^ 1 

Since the head of the expansion wave is at x ew = l/4\/2 for the logotrope, equation (|3.6|) 
shows that it reaches the surface of a collapsing cloud (of any mass) at 

tew = — — iff , (3.7) 

7T 

i.e., after ~1.80 initial free-fall times. At that point, equation ( |3.5D (or the ratio mo/m ew from 
eq. 2.22[ ) gives Mo/M tot — 0.03. If the similarity solution continues to hold after this time, then 
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the entire cloud will be accreted onto the core (i.e., Mo/M to t = 1 in eq. |3.5]) by 



2 /a/2\ 1/4 _ 

tend = ~ — t S ~ 4.32 iff . (3. 



7T \ TTlQ 



The inside-out collapse described by the expansion-wave solution is the gentlest possible, and 
therefore also the slowest. Solutions of the type ( |2.16| ), with C > \/2 for a sphere which is initially 
denser than allowed for equilibrium, lead to faster accretion because of an increase in mo (Table 
0) and because of an increase in the total mass (|3,3j), by a factor (C/\/2) 3 ^ 2 (see Appendix A). 



Realistically, once the expansion wave reaches the boundary of a cloud which is under a 
constant surface pressure, it should be reflected and driven back into the cloud as a compression 
wave. Thus, our similarity solution cannot be exact beyond t ew , and the resulting estimate of 
i cn d might be something of an upper limit. However, it does seem likely that more than one 
additional free-fall time is required to complete the collapse after i ew , because the accreting gas 



is in free is in free-fall only deep inside the expansion wave (at x < 0.1x ew ; cf. eq. p.l8f| ). In the 
absence of detailed numerical simulations, then, the most secure conclusion is that the collapse 
of a (near-)equilibrium logotrope lasts for t en d > t ew + ts ^ 3ffj after core formation. However, 
if the isothermal sphere is any guide, then the simulations of Boss & Black (1982) and Foster 
& Chevalier (1993) suggest that the correct time to Mq = Mt ot for a logotrope may in fact be 
fairly close to our equation Q3.8j ): The core in Boss & Black's numerical collapse of a non-rotating, 
equilibrium singular sphere grows at roughly the expansion- wave rate even after i ew ; and the 
results of our Appendix A can be used to compute a t enc j which is only ~20% longer than the total 
accretion time in Foster &; Chevalier's "£ max = 40" run for a cloud that is overdense by 10% and 
has high but finite concentration at t = —oo. 

Figure [3| shows the density profile and velocity field, in similarity space, for a full expansion- 
wave solution. The thick black lines lie inside the head of the expansion wave, which is at 
x ew ~ 0.177. At larger x, the density profile is given by the hydrostatic singular solution 
a = \f2jx. This is shown as the thin solid line in the top panel, where the continuous matching 
of densities across x ew is evident. The velocity field also joins continuously to the static v = at 
x ew . Table [3| lists a, v, and the mass m as functions of x inside the head of the expansion wave. 
The initial radius R of the cloud is fixed in real space (at least, until the expansion wave reaches 
it), and so over time moves steadily inwards in similarity space. The vertical lines in Fig. ^ are 
placed at X = R/att (eq. |3.6|| ), for the times t/ts shown. As discussed above, X = x ew after 1.80 
free-fall times. Finally, the dashed lines in both panels of the Figure are just the free-fall profiles 
of equation ( 2.17[) . Clearly, they do not describe the flow until the gravitational attraction of the 



core exceeds the pressure-gradient force, at x < 0.01 or so. 

A dimensional version of Fig. || is given in §4 below, but we now examine the dimensional 
core masses Mq and accretion rates Mq = 4(Mq) as functions of time in logotropes of various total 
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masses. These are obtained from equation (3.5): 
M = 5.41 x 10~ 6 M 



/M tot \ ( t s \ 4 / t 




4.8 x 10 5 yr I \10 5 yr 



\ -4 
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M = 2.16 x 10- iU M yr" 1 — ^ - - " — =— . (3.9) 

u J \M y ^4.8 x 10 5 yr y \ 10 5 yr/ y 1 

It will be recalled (from Table [l] above) that % = 4.8 x 10 5 yr is appropriate for a solar-mass, 
A = 0.2 logotrope with temperature T = 10 K and surface pressure P s = 1.3 x 10 5 k cm -3 K. 
The free-fall time of a logotrope with any other M tot , P s , or T can be computed as discussed in 
§2.1, so that equations ( |3.9| ) are quite generally applicable. The first of these relations is plotted 
in Fig. |] for a number of different initial cloud masses. In a C > \/2 collapse, our Mo and Mo 
should each be multiplied by (C/\/2) 3//2 (mo/6.67 x 10~ 4 ), with mo taken, for example, from Table 
H above. What was already apparent from equation ( |3.5[) is made more explicit here: a protostar 
in a collapsing logotrope gathers fully half of its mass Mo over the last 1 - (0.5) x / 4 ~ 16% of its 
accretion phase. By far the majority of a star's total formation time is spent waiting out a period 
of very slow and relatively unproductive infall in its parent cloud. 

Also shown in Fig. | is the line M = 0.155 M (t/10 5 yr)(T/10 K) 3 / 2 for a singular 
isothermal sphere of any total mass. There are at least three points of interest in a comparison 
between this and the logotropic curves. (1) Any isothermal sphere has the same time-invariant 
accretion rate, while more massive logotropes eventually accrete much more rapidly than clouds of 
lower M to t- This is because a softer-than-isothermal EOS generally gives rise to larger nonthermal 
motions on larger scales, i.e., to higher <j ave in higher-mass clouds (cf. §2.1). Roughly, then, since 
M ~ <7 3 /G, the mass originating from larger radii in bigger clouds will fall faster onto the core. 
(2) Generally speaking, low-mass stars take longer to accrete in a logotrope than in an isothermal 
sphere, while high-mass stars form on relatively shorter timescales. The extent of the discrepancy 
between the two models depends on the initial mass of the logotrope. Indeed, because of the effect 
noted in point (1), even the stellar mass for which the total accretion times coincide depends on 
Mtot- (3) This being so, there is a much smaller spread in star- formation times in logotropes than 
in isothermal spheres. This is due to the weak mass dependence in the free-fall time of a logotrope 
(Table [l]), and to the fact that M oc (t/t s ) 4 . Specifically, stars of mass O.3M < M < 3OM all 
tend to form within 1 — 3 x 10 6 yr, in logotropes of essentially any Mt Q t- In an isothermal collapse, 
where Mo oc t, the accretion times would differ by a factor of 100. 

It has to be recalled that these results are based on the assumption that the self-similar 
collapse solution continues to hold even after the expansion wave has been reflected from the 
boundary of a cloud. As discussed above, Fig. |I| may then represent upper limits to the true 
star-formation times in logotropes. On the other hand, this effect may not be very large; and in 
any case, it will be offset to some extent by the (outwards) radiation pressure from stars larger 
than a few M (Jijina & Adams 1996). Our analysis should then be a reasonably accurate 
first-order approximation to a rather complicated situation. 
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These points could have some bearing on the formation of bound clusters, which requires that 
a large fraction (^30%: Hills 1980; Lada, Margulis, & Dearborn 1984; Elmegreen & Clemens 1985) 
of the protocluster gas be converted to stars before the remainder is dispersed by supernovae or 
H II regions. In light of Fig. and the short accretion times of high-mass stars, this would seem 
to demand that high-mass stars start collapsing after the low-mass stars do — a conclusion which 
also has ramifications for the form and origin of the stellar IMF. 

Myers & Fuller (1992; also Myers & Fuller 1993; Caselli & Myers 1995) have introduced a 
"TNT" (thermal plus nonthermal) model of molecular clumps, which has an EOS that is explicitly 
isothermal at small radii, but becomes essentially polytropic (with 7 < 1) at large r. They use 
this model to estimate approximate accretion times for stars of various masses. Although our 
logotropic model is preferable to the TNT construction (because the latter has density singularities 
at the centers of even pre-collapse, t = —00 clouds), the Mq vs. t relation of equation ( |3.9| ) and 
Fig. |H may profitably be compared with that of Myers & Fuller. In particular, although our i enc j 
for clouds of various masses tend to be somewhat longer than the TNT estimates — which bypass 
an exact treatment of the expansion wave, and which are almost isothermal for low-mass clumps 
anyway — they have roughly the same order of magnitude. More importantly, the rather narrow 
range in accretion times, and the fairly rapid formation of high-mass stars, are common to both 
sets of results; indeed, these are generic features of any softer-than-isothermal EOS. 

4. Implications for Low-Mass Star Formation 

Our analysis has concentrated on spherical infall in a logotrope, and ignored the effects 
of mean magnetic fields (to whatever extent disordered fields contribute to the turbulence in 
interstellar clouds, they are already accounted for by the EOS), rotation and circumstellar disks, 
and radiation pressure from the forming star. All of these have been studied analytically or 
numerically for accretion in isothermal clouds (e.g., Galli & Shu 1993a, b; Tereby et al. 1984; Boss 
& Black 1982; Jijina & Adams 1996, and references therein). In addition, some connection has to 
be made with the issues of fragmentation and binary formation (see, e.g., Myhill & Kaula 1992). 
Clearly, the present study can only be viewed as one step towards a full theory of non-isothermal 
collapse. 

Nevertheless, some fairly robust conclusions can be drawn. In particular, although others 
have considered the effects of softer-than-isothermal, and even logotropic, equations of state on 
massive-star formation (e.g., Myers & Fuller 1992; Jijina & Adams 1996), the most important 
new point to be made here is that even low-mass (~ 1-M©) molecular clumps can be described by 
a logotropic EOS. Their collapse is therefore rather gentler than in the isothermal case, and their 
accretion rates must be much smaller at early times. 

Figure [5] shows the time evolution, from t = 5 x 10 4 to t = 8.6 x 10 5 yr, of the density profile 
and velocity field in an expansion wave moving through an initially hydrostatic A = 0.2 logotrope 
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with a total mass of IMq, a kinetic temperature T = 10 K (hence a c = 1.88 x 10 cm s _1 ), and 
a constant surface pressure P s = 1.3 x 10 5 k cm -3 K. The velocity scale at (eq. [ |2.10| ]) can be 
evaluated, at any time, given the surface pressure and the fact (cf. Table [l]) that P c /P s = 1.58 in 
the non-singular cloud at t = — oo; the dimensional r, p, and u of Fig. || then follow immediately 
from the x, a, and v of Fig. ||. At the point t = of core formation, the cloud is everywhere at rest 
and follows the singular density profile p(r) = 1.3 x 10 -20 g cm -3 (r/i?) -1 . At later times, this 
profile still holds outside of the expansion wave; it is sketched as the broken line in the top panel 
of Fig. ||. For reference, we also have that the initial free-fall time of this cloud is tg = 4.8 x 10 5 
yr (Table |IJ); the initial ratio of central and mean densities is p c j p ave = 4.20; and the initial radius 
is R = 2.9 x 10 17 cm (eq. |3.4j| ; this is within a factor 2 of the radius of a solar-mass isothermal 
sphere under the same surface pressure). We have only plotted the densities and velocities for 
t < 8.6 x 10 5 yr = 1.80 tg, i.e., prior to the time t ew when the head of the expansion wave meets 
the initial boundary of the cloud. 

A cloud which is initially rotating at a constant rate f2 will depart from the spherical 
expansion- wave solution at small r <^ Rc, where Rc is the centrifugal radius at which the circular 
speed (GMo/r) 1 / 2 just equals the infall velocity —u. For an M tot = 1M@ logotrope with the 
surface pressure and temperature we have assigned here, it is given by 

n 2 M in / n \ 2 / t \ 4 , . 

Rc = — = 3.0 x 10 10 cm hr-TT-T T7^- • (4.1) 
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The first of these equalities comes from Jijina & Adams (1996; note the slightly different notation); 
the second, from our equation ( |3.9| ), along with A = 0.2 and the fact that P c = 1.58P S for this 
specific example. Thus, a disk of 100 AU in extent will be built up within ~1.5 x 10 6 yr in our 
model, as compared to ~7.5 x 10 5 yr for a singular isothermal sphere at T = 10 K. Scaling of this 
result for clouds with different total masses, surface pressures, and temperatures involves a change 
in the coefficient of Mq, as described after equation (|3.9|) , and a change in P c /P s from, e.g., Table 
|]. In any case, the magnitude of Rc at a specified age t gives some indication of the radial range 
in which the results of Fig. |E| are valid in an initially rotating cloud. Calculation of the analogous 
length scale for a magnetized sphere — the so-called magnetic radius of Galli & Shu (1993a, b) - 
will also be important, but is beyond the scope of this paper. 

There are two main points to Fig. [El. The first is that the densities far behind the head of the 
expansion wave are well in excess of the initial hydrostatic values (see also Fig. |3| above), and are 
increasing with time. This is in stark contrast to an isothermal collapse, where p{r) falls below its 
equilibrium value and decreases with time at small radii (cf. Fig. 3 of Shu 1977). The difference 
stems in part from the increasing Mq in our case, but more fundamentally from the fact that the 
initial r~ l density profile of the logotrope is shallower than the final r -3 / 2 which must obtain near 
the central core; there is relatively more mass in the outer regions of the static cloud than in the 
collapsing one, so a backlog of sorts develops at small radii. The opposite is true of a singular 
isothermal sphere, which begins with a static r -2 density profile. More quantitatively, the scaling 
a oc x -3//2 implies p(r, t) oc r _3//2 t for the logotrope, but p(r, t) oc r^ 3 / 2 t -1 / 2 for an isothermal 



-18- 



EOS. 

This situation calls to mind observations of the far-infrared emission of the NGC 2071 region 
in Orion B (Butner et al. 1990). Radiative-transport models are consistent with a central star of 
around 5M there, and a p oc r~ 3 / 2 density profile at radii r > 2000 AU. However, the value of 
the density at 2000 AU given by Butner et al.'s radiative transport analysis (and also suggested 
by independent molecular data) is at least one full order of magnitude larger than that expected 
for an isothermal expansion-wave collapse. One explanation for this is that 2000 AU is deep inside 
the expansion wave for an isothermal collapse age of ~ 5 x 10 5 yr (which follows from an assumed 
T = 35 K and Mq = 5Mq), such that significant depletion would have occurred there. There is 
no such depletion, but rather an enhancement, in our logotropic model, and these data could 
therefore provide an important test of it. However, a proper analysis of NGC 2071 would require a 
self-consistent treatment of radiation pressure and rotation, as well as a reliable feel for the initial 
total mass involved, and is too detailed to be undertaken here. 

The second point of Fig. |5| is that, naturally, the infall velocity at any given radius also grows 
over time, but at a rate which is quite different from that in an isothermal sphere. In particular, 
since it oc (Mo/r) 1 / 2 at small radii, we have that u oc r~ l l 2 t 2 in a logotrope, while u cx r" 1 / 2 ^/ 2 
in the isothermal case. At any fixed r, the velocities in the two models will agree at some t; the 
logotropic u must then be smaller than the isothermal one at earlier times, and larger later on. 
With Mtot = lAf© specifically, a comparison of Fig. | and Shu's (1977) Fig. 3 shows that we are 
firmly in the former, lower- velocity regime at all radii and throughout the entire duration t < t ew . 

For example, let us focus on the scale r = 10 16 cm ~ 0.003 pc. Our model gives infall 
velocities there (and at larger radii) which are less than the thermal linewidth, a c ~ 0.2 km s , 
for all t < 8.6 x 10 5 yr. By contrast, a solar-mass, singular isothermal sphere attains —u ~ a c at 
the same radius by t = 6 x 10 4 yr. This is one example of how the systematic velocity field at 
observable radii in a logotrope of any mass will be effectively swamped by random velocities for a 
significant fraction of the duration of its collapse; even a molecular clump which is well into the 
accretion phase can appear to be in virial equilibrium. Because the accretion rate increases with 
time, at any instant the majority of clumps will be in a fairly inconspicuous evolutionary state. 
This could explain the apparent paucity of obviously collapsing regions of the ISM. Similarly, for 
those clouds which do show observable signatures of collapse, the adoption of a logotropic EOS 
leads to larger inferred dynamical ages than an isothermal model would suggest. 

Currently the best direct, kinematic evidence for collapse is in the isolated Bok globule B335, 
where molecular line profiles appear consistent with infall velocities — u ~ 1.2 km s -1 at radii 
r ~ 5 x 10 15 cm (Zhou et al. 1993). Given that — u = (2GMo/r) 1 / 2 deep inside any expansion wave 
(which radial scaling is also consistent with the B335 data), this implies the existence of a central 
protostar (or star plus disk) with mass Mq ~ 0.25M Q . As well, 1.3-mm continuum measurements 
of dust emission indicate a circumstellar envelope of about O.8M0 (Bontemps et al. 1996). 
The physical parameters which we have specified in deriving Fig. || therefore seem particularly 
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suited to a discussion of B335.Q In particular, model-independent considerations suggest that 
Mo/M to t ~ 0.25 in this object. If it is undergoing a logotropic collapse, it must then be older than 
i ew = 8.6 x 10 5 yr, because our similarity solution gives a much smaller M /M tot = 0.03 at that 
time. (In addition, Fig. || shows that — u is rather less than the observed value at r = 5 x 10 15 
cm for all t < t cw .) That is to say, the expansion wave should already have been reflected from the 
boundary of B335. 

If B335 is much older than t ew , then it will no longer follow its original r~ l density profile; the 
entire cloud should be taken up in the collapse and the free-fall p oc r~ 3//2 should hold essentially 
throughout. In fact, this is consistent with the study of Zhou et al. (1990), who conclude that the 
densities in the inner regions of the globule are best described by an r~ L5 power law. (As a test 
of the isothermal collapse theory specifically, they go on to argue that at large radii the densities 
are consistent with p oc r~ 2 . However, there is no direct evidence for this, or for any other power 
law but —1.5; the data are "insufficient to determine independently the parameters of a density 
distribution with two power-law regimes" [Zhou et al. 1990].) 

Just how old is B335? Equation (|3.9j ), given Mq = 0.25M©, predicts a collapse age of 
t = 1.5 x 10 6 yr (if our fiducial surface pressure is correct here), and hence a time-averaged 
accretion rate of (Mo) = Mo/t = 1.7 x lCP 7 A/0 yr -1 . These numbers are, moreover, quite robust 
again any uncertainties in the initial total mass: even for an unrealistically extreme M tot = WMq 
(in which case Mq/M^ ~ 0.03 and the expansion wave would have just reached the edge of the 
cloud), we find t ~ 1.2 x 10 6 yr. By contrast, Zhou et al. (1993) discuss their data in terms of 
the isothermal expansion wave model and estimate t = 1.5 x 10 5 yr, Mq = 2.8 x 10 -6 Mq yr -1 , 
and Mq = OAMq. It seems inevitable that our model requires an age for B335 which is nearly an 
order of magnitude older, with a time-averaged accretion rate that much smaller, than previously 
thought. 

B335 is an example of the low-mass "Class 0" sources recently identified by Andre et al. (1993) 
and Andre & Montmerle (1994). These typically have (cold) single-blackbody spectral energy 
distributions, but are invisible at short wavelengths A ^ 10 pm. This is attributed to veiling by 
a circumstellar envelope whose mass exceeds that of the central protostar — a conclusion which 
is supported by direct, mm-continuum measurements of the envelope masses — and therefore 
implies that these objects are in the earliest stages of collapse yet to be observed. Indeed, the 
theory of isothermal accretion applied to the Class sources generally implies ages of a few 10 
yr (e.g., Barsony 1994; Pudritz et al. 1996). Although B335 might then be somewhat older than 
average even in the context of isothermal models, it seems likely that the relative age increase 



2 In fact, Zhou et al. (1990) favor a kinetic temperature of 13 K for B335. A T — 13 K logotrope with no mean 
magnetic field has a critical mass of 155M© and a free-fall time of 6.6 x 10 5 yr, by equations ( |2.2j ) and (2.3). 
Thus, B335 has M to t/M C rit — 0.0064, and hence a truncation radius R/ro — 1 from Table |l[ This in turn implies 
iff /iff.crit — 0.74, or tg = 4.9 x 10 5 yr. Our timescale estimates, which strictly apply to a T = 10 K solar-mass 
logotrope with tg = 4.8 x 10 5 yr, are then relevant to B335 as well. 
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we deduce for it should also be applied to other Class sources. Note, though, that most of this 
age increase is required to accomodate an early period of very gentle infall, with small central 
source masses and low accretion luminosities. Most of a 1M@ star is accreted over a period of only 
~ 2 — 3 x 10 5 yr in a logotrope. This is much smaller than the total collapse time (cf. Fig. |j) and 
is still consistent, for example, with empirical estimates of the formation time of visible T Tauri 
stars in the Taurus- Auriga complex (Myers & Benson 1983; Myers et al. 1987). 

Finally, if our low mean accretion rate in B335 is typical, then it should have some bearing on 
the "luminosity problem" of young stellar objects in Taurus- Auriga. Very briefly, the bolometric 
luminosities of Class I sources there (which are embedded in less massive circumstellar envelopes 
than Class 0's, and hence appear to be in later stages of collapse; Andre &: Montmerle 1994) 
are generally about an order of magnitude smaller than expected from number-count arguments 
(Kenyon et al. 1990; Kenyon et al. 1994) and fits to individual spectral energy distributions 
(Kenyon, Calvet, & Hartmann 1993; see also Adams, Lada, & Shu 1987). A similarly low 
luminosity has also been noted for the prototypical Class source VLA 1623 (Pudritz et al. 1996). 
The problem derives, in large part, from the isothermal assumption of a time-independent 
accretion rate Mq, and Kenyon et al. note that it could be alleviated if stars accumulate most of 
their mass on timescales which are short compared to the total duration of their embedded phase. 
We have seen that such a situation arises naturally in the collapse of a logotrope. A revision of the 
assumed EOS in molecular clouds may therefore play an important part in the eventual resolution 
of this issue. 

5. Summary 

McLaughlin & Pudritz (1996) have shown that a logotropic equation of state, 
P/P c = 1 + A ln(p/p c ) and A = 0.2, provides the most self-consistent description of the 
internal, equilibrium structure of both low- and high-mass dense clumps within GMCs, and of the 
global properties (Larson's laws) of the giant cloud complexes themselves. We have now used this 
model to investigate the collapse of interstellar clouds. There are a number of consequences for 
issues of not only high-mass, but also low-mass, star formation. 

We have derived similarity solutions for the collapse of a gaseous logotrope. These solutions 
are analogous to those of Shu (1977) for the collapse of singular isothermal spheres; in particular, 
an expansion-wave solution — the inside-out collapse of a cloud which is initially in virial 
equilibrium — is identified as the one of most physical interest. The predicted density profile is 
p oc r _1 for gas which has not yet been reached by the expansion wave, and p oc r~ 3 / 2 far behind 
the collapse front, where infall velocities approach their free-fall values — u = (2GMo/r) 1//2 . These 
latter scalings are model-independent, as they follow directly from the presence of a central point 
mass (protostar) Mq. However, the core accretion rate in a logotrope increases with time (as 
Mq oc i 3 , so that Mq oc t ), while it is constant for an isothermal sphere. 
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This implies that Mq in a logotrope will be smaller than the isothermal value at early times, 
and larger at late times. Low-mass stars therefore take longer to form, while high-mass stars are 
built up much more rapidly. With % the mean free-fall time of a logotrope (typically of order 
5 x 10 5 yr), our similarity solutions are strictly valid for t < 1.8%. At that point the expansion 
wave reaches the edge of a cloud — with only 3% of the total initial mass having been accreted 
onto the core — and hydrodynamical simulations will likely be required to exactly describe the 
later stages of evolution in which most of the final mass of a star is accreted. Nevertheless, a 
simple extrapolation of our present results suggests that if a cloud of any mass were to completely 
collapse to a single star, it would do so by t = 4.3 %. A weak dependence of tg on mass in 
logotropes, together with the result Mq oc i 4 , guarantees that the accretion times of stars of widely 
disparate masses differ only slightly: stars of anywhere from 0.3 to 3OM are all expected to form 
within ~ 1 — 3 x 10 6 yr. This conclusion, which is similar to that reached in other studies of 
non-isothermal collapse, has important implications for cluster formation and the stellar IMF. 

Another clear difference between the logotropic and isothermal models lies in the evolution of 
the density at small radii in a cloud: p ultimately increases with time at any given radius in our 
model, but would decrease in a singular isothermal sphere. Infall velocities u, on the other hand, 
grow with time in any collapse model, but the rates at which they do so are quite dependent on 
the assumed equation of state. Early on in its collapse, the logotrope can have much smaller u 
than an isothermal sphere. It therefore takes longer for observable systematic velocities to develop 
at observable radii within a cloud, and collapse will initially be harder to detect directly. This 
implies that (1) early on in its accretion phase, a logotrope will effectively appear to still be 
in virial equilibrium, such that the number of molecular clumps which are in fact undergoing 
collapse could be much larger than previously suspected; and (2) current estimates for the ages of 
(low-mass) young stellar objects, which are based on isothermal collapse models, should be subject 
to a significant upwards revision. In the Bok globule B335 specifically, we find an age which is 
certainly greater than 8.6 x 10 5 yr, and probably ~1.5 x 10 6 yr, older by perhaps an order of 
magnitude than suggested by the isothermal expansion-wave theory. This result could have wider 
implications for the ages of Class protostars (of which B335 is one) in general. 

This work was supported in part by the Natural Sciences and Engineering Research Council 
of Canada. 

A. COLLAPSING POLYTROPES 

A.l. Similarity Solutions 

In many ways, the isothermal sphere and the logotrope are profitably viewed as opposite 
extremes of the family of negative-index polytropes, which have P oc p 1 with < 7 < 1. Although 
incomplete as models of real interstellar clouds, these have been considered, for example, by 
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Maloney (1988), McKee & Zweibel (1995), and MP96. The case 7 = 1 gives an isothermal 
equation of state, while most features of the logotrope are recovered by formally setting 7 = 
in what follows. The "negative-index" label refers to standard polytropic notation, which has 
7 = 1 + 1/N. An isothermal EOS then corresponds to N —* —00; the logotrope, to N —* —1. 



To derive similarity solutions for collapsing polytropes, we write the EOS as 

P = Kpi . 



(Al) 



The similarity variable x, the dimensionless density a, the velocity v, and the mass m are all 
defined as in equations ( plf ) and ( |2.9| ) above; and since dP/dp = Kjp' 7 ^ 1 , the appropriate velocity 
scale is just 



at 



1/2 



(A2) 



As discussed in §2.1, we are concerned with the evolution of a cloud only after the time t = 
when a core has formed at its center. The density structure is then singular throughout, and K is 
a constant in space and time, with the value cr^pj -1 it had in the initial (t = —00) cloud. 



If we write 



then equations (2.5) 



a 



7-1 



P = (2- 
) give the following: 
d lna 
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.7-1 



dx 
dv 
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P + (1 - 7) W + v] 
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P 



2a 1 - 



(3 + (1 - 7) [2a 7_1 + v/3 



ax 
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(A3) 

(A4) 
(A5) 

4-37- (A6) 

Setting 7 = 1 in expressions (A2) onwards recovers the results of Shu (1977) for the isothermal 
sphere; if instead 7 = (but is replaced by AP C ), then the equations of our §2.2 are obtained. 
No constraint has yet been imposed on the value of 7, and in particular it could well be > 1 
at this point. In this regime, treatment of the case 7 = 4/3 evidently requires some care. We 
do not discuss this point further, but note that it relates to the instability of positive-index 
polytropes with 1 < 7 < 4/3 (e.g., Chandrasekhar 1967). Goldreich & Weber (1980) have studied 
the homologous collapse of a 7 = 4/3 polytrope, in connection with the evolution of degenerate 
iron cores in supernovae. More generally, Suto & Silk (1988) solve equations for the self-similar 
collapse of polytropes with arbitrary indices 7. The physical requirement K = constant identifies 
our equations (A4) and ( |A5| ) with equations (15) of Suto & Silk, in their "n = 2 — 7" case. 



Still without restricting the value of 7, equations (A4) - ( |A6| ) admit two analytic solutions. 
One of these is the same uniform-density case given by equation Q2.14| ); the other is the static 
solution, 

2 ~T ^ ,(4-3r)/(2-7) 



a = Chse x 



-2/(2- 7 ) 
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4-37 



Chse x ( 



(A7) 
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where the density normalization is 



Chse 



2(4 - 37) 
(2- 7 )2 



1/(2-7) 



(A8) 



The singular density profile thus obtained is an exact solution of the equation of hydrostatic 
equilibrium for a self-gravitating polytrope of any (negative) index (cf. MP96). Again, if 7 = 
here, the logotropic result ( |2.15 ) obtains, while if 7 = 1 we have a = 2/x 2 and m = 2x, in 
agreement with Shu (1977) for the singular isothermal sphere. 

In the context of interstellar clouds, we are interested only in the range < 7 < 1 — and 
particularly the 7 = limit of the logotrope — since the velocity dispersion P/p then increases 
with radius. For such 7, equations (A4) - ( |A6| ) in the limit x — > (late times or small radii) and 
\v\ S> x give 

0) nwmo ' (A9) 

which are minus solutions. In the opposite limit x — > 00 (early times or large radii) and \v\ <C x, 
we have instead 



a 



(4 - 3 7 ) 



a 



Cx- 2 /(2-7) 



•7 



4-37 



c 



x 



-7/(2-7) 



m 



2 ~ 7 Cx (4-3 7 )/(2- 7 ) 
L-37 

(A10) 

The density normalization C is essentially a free parameter, but it must exceed the hydrostatic 



equilibrium value Chse (eq. | A8|) for inflow (v < 0). As usual, then, this family of solutions 
represents the collapse of non-equilibrium spheres which are everywhere denser (by a factor 
C/Chse) than a virialized one with the same volume. There is also a one-to-one relationship 
between C and mo for a complete solution which extends from x = 00 (the moment of core 
formation) to x = (the end of accretion), and the limit C — > C^ SE corresponds to the 
expansion- wave solution. (Note that these asymptotic solutions differ from those cited by Larson 
[1969], which hold for plus solutions where \v\ increases at large x but is still <C x. A sufficient, 
though not necessary, condition for the existence of such solutions is 7 > 1/3.) 



A. 2. Accretion Rates 

The reduced core mass mo translates, by equations (|A2j) and ( p.9[ ), to a dimensional 

where tg is the free-fall time evaluated at the initial mean density p a ve; and where we have used 
the fact that K = a^p^T 1 . By construction, the linewidth a c at the center of a cloud is taken to be 
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purely thermal — a 2 = kT/^niH — and the EOS automatically accounts for nonthermal motions 
elsewhere. In any case, the mass accretion rate is therefore 

3/2 /3^^\ 3(1 - 7)/2 m 3(H W 



M = (4-3 7 ) 7 ^^--^-j ^ = (4-3 7 )(M ), (A12) 

and only the singular isothermal sphere has an accretion rate Mq = (Mo) = mocr^/G that is 
constant in time. Formally, the results of §3 above are recovered by substituting A 3 / 2 for 7 3 / 2 and 
setting 7 = elsewhere in these expressions. 

Now, by manipulating the results given in Appendix C of MP96, the total mass and radius of 
a singular pressure-truncated polytrope can be written as 

3 7 \ 3/2 /4_3 7 \ 37/2 / C x 3/2 . x 3(l- 7 )/2 n 3 



f '£3 Us it) ^ <-> 



and 

^f) (s^) (as) (£) (A14) 



7T 



which should be compared to equations ( p.3| ), (3^4) in the logotropic limit 7 = 0. The factor 



C/Chse here is required for consideration of the collapse ( A10| ) of a sphere which is everywhere 
denser than a virial-equilibrium configuration with the same radius R and temperature a c . Where 
the expansion wave is concerned, of course, the initial cloud is taken to be exactly virialized, and 
C = Crse- 

Given these results, the central point mass builds up as 

Mo 7T /V 4- 3 7 \ 3(1_7)/2 / 1 C V 3 / 2 /t\ 4 ~ 37 , A1C . 

m , (A15) 



Aftot 8\82- 7 ; V2- 7 C HS e/ Vfe 

and the boundary X = R/att is given by 

v 4/8 2- 7 V 1 ~ 7)/ V i c \ 1/2 fty- 2 , A . 



If our similarity solution is taken to hold at all times, then equation ( |A15 ) gives the time to 



accrete an entire cloud onto the core (Mo/M to t = 1)- For the expansion wave, with C = Chse, 
this ranges between t ent ± ~ 2.61 tg for 7 = 1, and t cnc i ~ 4.32 tg for 7 = 0. (The relation between 
7 and mo is given by Table |I] in Appendix B.) Strictly speaking, however, boundary effects are 
important at late times in finite spheres; in particular, the expansion wave will reach the edge of a 



cloud when X = x ew (eqs. P3| , |B4|), i.e., at 



K/ - 7 t s < W. (A17) 
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B. CRITICAL POINTS 



Collapse flows in polytropes with 7 < 1 (where again 7 is set to for the logotrope) may 
involve passage through a critical point, x* > 0, where the velocity v* < and density a* satisfy 

„(7-l)/2 



Equations 



(2 - 7 )x* 
and (|A!3) then demand that 



a* 



+ (1 -i)a 



(7-l)/2 



-a 



7-1 



4 — 37 

in order for (da/dx) X:t and (dv/dx) Xt to be finite. 

The critical point x ew for which v* = has, by ( [Bl| ) 



(l-7)(2-7)x, 



(Bl) 



(B2) 



a h - 1)/2 



so that (|B2j ) gives 



a ew = 2(4 - 3 7 ) 



2(4 - 3 7 ) 



1 V(2-7) 



-2/(2-7) 



(B3) 



(B4) 



(2 - 7 ) 2 

which agrees with the static solution (eqs. [|A7f| , |A8|| ) of the fluid equations at x ew . The collapse 
solution which begins with v — ► — 00 at x = 0, and goes on to pass through x cw with v = 0, may 
therefore be continued to larger x by attaching the static solution. This composite flow is just the 
expansion wave. 

We now confine our attention to critical points x* < x cw , so that v* < and a* > a ew - In 
general, given x*, it is a simple matter to find a* and v* from equations (]B^ ) and (Bl). We then 
set 

a(x) = a* + ai(x — x*) + 0[(x — x*) 2 ] and v(x) = + 61 (x — x*) + C[(x — x*) 2 ] 

in some neighborhood of x*. Substitution of this first-order expansion in equations (A4) and (A5) 
leads to a pair of coupled quadratic equations for the slopes a\ and b\. Under the constraint (B2), 
these yield a single cubic equation for a±, and thus three eigensolutions for the flow through a 
critical point. The first of these satisfies the simple relation 



61 = (2 - 7) + 



1 



(B5) 



which serves only to move the flow along the locus ( |B1| ) of critical points. Such a propagation of 
singularities is so contrived as to be unphysical, and we therefore discard this solution. This leaves 
us with 
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where 

fei - -(9- 77) + -a^ 

fc 2 ee a, + 2(l- 7 )(5-3 7 )-fcMal^ 1)/2 + 4«r 1 - (B7) 

For < x ew and < 7 < 1, ai is always real and negative. If the '— ' sign is taken to 
solve for it in equation (|B6|), we find b\ > 0, and such solutions have v more negative at smaller 
x. These therefore correspond to the minus solutions of Shu (1977), or the type 2 solutions of 
Hunter (1977) and Whitworth & Summers (1985). In the isothermal limit of 7 = 1, the result is 
a\ = —2/xl and b\ = 1/x*. If instead the '+' sign is taken in equation (|B6[), then 61 < and we 
recover the plus (or type 1) solutions; in the isothermal case, a\ = 2/x* — and b\ = 1 — 1/x*. 
The plus solutions for this or any other 7 have \v\ growing for x > x*, and vanishing at some 
nonzero x < x*. In addition, they are underdense at large x, relative to the singular hydrostatic 
profile, and may therefore be of interest as (time-reversed) wind solutions. 

The slopes a\ and b\ at be found analytically, given any 7, for both the minus and 

the plus solutions which flow through that point; again, the minus solution at x ew is the expansion 



wave. Using equations (B3) and (B4) in (|B6|) and (|B7|), we have the following: 



• 3/5 < 7 < 1, minus solution: 

= 37-5 a (3_ 7)/2 < Q ^ 6i = 57-_3 >(K 
1+7 1+7 

Although CK ew and v ew match the density and velocity of the corresponding static solution at 
s ew , the derivatives do not agree there (in particular, dv/dx = for the static outer cloud). 
The density profile and velocity field of an expansion wave in this range of 7 — which 
includes the isothermal sphere — are therefore discontinuous at x cw - 

• 0<7<3/5, minus solution: 

a x = -2af- l)/2 < and 61 = . (B9) 

Now these (left) derivatives do match the (right) derivatives of the static solution at the 
head x ew of the expansion wave (see eq. [|A7|), and the full flow is everywhere continuous. 
The logotrope corresponds to 7 = 0, and therefore shows this behavior. 

The plus solutions at x ew for 7 > 3/5 have the same a\ and 61 as the minus solutions for 7 < 3/5; 
and the plus solutions at x ew for 7 < 3/5 have the same derivatives as the minus solutions for 
7 > 3/5. 



Equations (|B8|) and (B9) also show that, for 7 < 1, the velocity slope b\ of the expansion 



wave at x ew is always shallower than the gradient dv*/dx* of the locus (|B1|) of critical points there, 
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and the infall velocity of the material just inside the head of the wave is less than the critical 
value. However, \v\ increases without bound towards x = for a minus solution, so any gas EOS 
which is softer than isothermal (including logotropic; §2.3) leads to an expansion-wave flow that 
eventually crosses a second critical point x* < x ew . Since the expansion wave also stands as the 
limit of out-of-equilibrium collapse solutions (C > Chse; Appendix A), some of these must also 
pass through two critical points. By contrast, in the isothermal case the expansion wave is just 
tangent to the locus of critical points at x ew = 1 ; and all C > 2 infall solutions avoid it altogether. 

Finally, it is clear that the head x ew of the expansion wave; its second, inner critical point x*; 
and the core mass/central accretion rate tuq are all uniquely defined by the specification of the 
polytropic index 7. These connections are shown in Table |j. Also given there is the total mass 
m ew interior to x ew (from eq. |A7| ). The ratio mo/m ew is the fraction of infalling material which 
is already taken up in the central core itself. 
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Table 1. Properties of Pressure- Truncated A = 0.2 Logotropes. 



R/r 
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7.2 
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0.0149 


9.3 


87 


6.7 


24.37 b 


115.58 


20.00 


0.0130 


9.0 


92 


5.8 



a Calculated assuming no mean magnetic field; surface pressure P s = 
1.3 x 10 5 k cm -3 K; and kinetic temperature T = 10 K. 

br Truncation radius of a critically stable cloud, M tot = M cr i t . 
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Table 2. Core Mass vs. Density Normalization. 
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Table 3. Logotrope Expansion Wave. 
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Table 4. Polytrope Expansion Waves. 
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Fig. 1. — Density profiles of a pressure-truncated logotrope and the Bonnor-Ebert isothermal 
sphere. Both configurations are in hydrostatic equilibrium, and correspond to critical-mass clouds, 
i.e., an increased central concentration in either case would destabilize the equilibrium. The scale 
radius is r$ = 3cr c /(47rG / o c ) 1//2 , and p s is the density at the cloud surface. The broken lines trace 
the singular profiles for each equation of state; note how such a solution mirrors the true behavior 
of the logotrope at most r, but bears little resemblance to any stable isothermal sphere. 

Fig. 2. — Velocity fields for various logotrope collapse solutions. Solid curves are minus solutions, 
and broken lines are plus solutions; the dash-dot line traces the locus of critical points. The heavy 
solid line represents the expansion wave, which continues to x = oo with v = and densities given 
by the hydrostatic singular solution. The minus solutions to the left of the expansion wave reach 
v = at some finite x, and therefore have vanishingly small spatial extent at t = 0. The minus 
solutions to the right of the expansion wave follow the collapse of spheres which are initially out of 
hydrostatic equilibrium. The expansion wave is the limit of both sequences. 

Fig. 3. — Density profile and velocity field, in similarity variables, of the expansion-wave solution 
for a collapsing logotrope. The heavy black lines trace out the region of dynamical collapse inside 
the expansion wave, like the same line type in Fig. ^, while the thin solid line in the top panel is 
the hydrostatic singular solution a = 2 1 / 2 /x. The dashed lines are the expected free-fall solutions 
a ~ £~ 3 / 2 and v ~ x -1 / 2 . Vertical (dotted) lines show the time-dependent location X = R/att of 
the cloud's outer radius. 

Fig. 4. — Expansion-wave accretion timescales for stars in an A = 0.2 logotrope and a singular, 
T = 10 K isothermal sphere. Low-mass stars take longer to form in the logotrope, while more 
massive ones accrete more rapidly. 

Fig. 5. — Time evolution (from left to right), in dimensional variables, of the density profile and 
velocity field for a collapsing logotrope of 1M© with a kinetic temperature T = 10 K and under 
a surface pressure P s = 1.3 x 10 5 k cm -3 K. The solid curves lie in the interior of the expansion 
wave (and therefore correspond to the heavy black lines in Figs. [2] and |3|), and the broken line in 
the top panel traces the hydrostatic r _1 density profile out to the boundary R of the initial cloud. 
The times shown are roughly the same t/tg used to locate X in Fig. [3|. 
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